First and second-type self-similar solutions of implosions and 
explosions containing ultra-relativistic shocks 
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ABSTRACT 

We derive self similar solutions for ultra-relativistic shock waves propagating 
into cold material of powerlaw density profile in radius p oc r~ k . We treat both 
implosions and explosions in three geometries: planar, cylindrical and spheri- 
cal. For spherical explosions these are the first type solutions of Blandford and 
McKee for k < 4 and the second type solutions found by Best and Sari for 
k > 5 — a/3/4. In addition we find new, hollow (with evacuated interior), first 
type solutions that may be applicable for 4 < k < 17/4. This "sequence" with 
increasing k of first type solutions, hollow first type solutions, and then second 
type solutions is reminiscent of the non-relativistic sequence. However, while in 
the non relativistic case there is a range of k which corresponds to a "gap" - a 
range in k with neither first nor second type solutions which separates the hollow 
first type solutions and the second type solutions, here there is an "overlap": a 
range of k for which current considerations allow for both hollow first and second 
type solutions. Further understanding is needed to determine which of the two 
solutions apply in this overlap regime. We provide similar exploration for the 
other geometries and for imploding configurations. Interestingly, we find a gap 
for imploding spherical shocks and exploding planar shocks and an overlap for 
imploding planar solutions. Cylindrical configurations have no hollow solutions 
and exhibit direct transition from first type to second type solutions, without a 
gap or an overlap region. 
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1. Introduction 

Self similar solutions provide some of the greatest simplifications to one dimensional 
flows. Self similarity allows the reduction of the partial differential equations, which contain 
two independent variables (space and time), into a set of ordinary differential equations, 
where the single independent variable is a combination of space and time. The ODEs are 
then relatively easy to solve numerically or even analytically in some cases. They describe 
the asymptotic behavior of one dimensional flows in a variety of circumstances, typically 
far away from the initial conditions and provided that the boundary conditions contain no 
spatial scale (Some exceptions apply. For example, self similarity can prevail in exponential 
density gradient in planar geometry). 

Perhaps the most famous hydrodynamic self similar solution is that due to Sedov(l), 
Von Neumann(2) and Taylor(3) known as the Sedov- Taylor solution. It describes an 
explosion in which a strong shock wave propagates into cold surroundings whose density 
profile decreases as 1 p oc r~ k . They used conservation of energy to obtain the scaling of the 
shock radius as a function of time. Such solutions are called first type solutions. Zeldovich 
and Raizer (4) discuss a self similar solution describing imploding shock waves in a constant 
density environment. In contrast to the strong explosion problem, energy considerations 
cannot be used to deduce the scaling of the shock radius as a function of time. Instead the 
scaling of the radius as function of time must be found by demanding that the solution pass 
through a singular point of the equation. Such solutions are called self similar solutions of 
the second type. More recently, Waxman and Shvarts(5) showed that if the density falls 
fast enough (k > 3), energy considerations give the wrong scaling. They also showed that 

Tn the non relativistic literature, the density powerlaw index is usually denoted by uj 
while in the relativistic temperature the letter k is used. Here we use the letter k for both 
as our focus is the relativistic regime 



-4- 



the solution should be of the second type for k > 3.26. 

Recent astrophysical discoveries, mostly the afterglows of gamma ray bursts, led to 
an increased interest in the relativistic analogs of these solutions. We find that there is a 
considerable similarity between the relativistic and the the non relativistic regimes. The 
relativistic version of the Sedov- Taylor solution was found by Blandford and McKee in 1976 
[BM hereafter, (6)]. They provide the solution to the ultra relativistic strong explosion 
problem in spherical geometry for density profiles with k < 4. Best and Sari(7) have found 
solutions of the second type for k > 5 — a/3/4. Perna and Vietry(8) found relativistic 
solutions for shock waves propagating in an exponential density gradient. These are the 
relativistic analogs of the non relativistic solutions first found by Raizer(9). 

In this paper, we provide a more complete exploration of the possible self similar 
solutions in the ultra-relativistic case. We provide solutions in three geometries, planar, 
cylindrical and spherical for both imploding and exploding shock waves. Our investigation 
reveals some unexpected puzzles. For example, for spherical explosions, we find a range 
of density profiles 5 — a/3/4 < k < 17/4 for which current considerations allow for both 
first and second type solutions. Understanding which of these solutions will apply is left for 
later research. Fortunately, all our solutions are analytic. This is an advantage over the non 
relativistic case where the self similar solutions are either implicit (first type) or numerical 
(second type). This property may facilitate the understanding of the interplay between first 
and second type solutions. 

The plan of this paper is as follows. We begin in section §2 by writing down the one 
dimensional flow equations in planar, cylindrical and spherical symmetries. We then take 
the ultra relativistic limit and assume self similarity so the partial differential equations are 
reduced into ordinary differential equations. In section §3 we explore the range in density 
powerlaw index for which first type solutions, which obey global conservation law, apply 
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and contain a finite energy. In section §4 we find second type solutions and explore the 
regimes where those are valid. 



2. The self similar equations in arbitrary dimensions 

The flow equations are given by the conservation of energy, momentum, and particle 
number. Since in the relativistic case the particle number density does not play a role in 
the dynamics of the system, the conservation of particles equation is decoupled from the 
other two equations. Energy and momentum conservation yield 

^ + P 2 P) + ^l^[3(e + P ) = (1) 
^(e + p) + ^A¥(e + p) + ^ = 0, (2) 

and the particle conservation equation is 

— in +— — r Q 7/ 3n = 0, (3) 
dt 1 r a dr w 

where a — 0, 1, 2 for planar, cylindrical and spherical symmetries respectively. 

We now assume that the flow has a characteristic Lorentz factor T and a related 
characteristic position R. R evolves with time in a way that describes motion with a 
Lorentz factor T. In the extreme relativistic case this implies R = - l/Y 2 « 1 - l/2r 2 . 
For the problem of interest here, it is natural to choose the position of the shock, and its 
Lorentz factor as R and V. The characteristic length scale of the flow behind the shock 
must be R/T 2 , so the similarity variable is 

C = ^ = ni-r/R). (4) 

BM doubted the use of equation (4) for k > 4; however, Best and Sari (7) demonstrated 
that this definition of the self similar variable is applicable in general. 



- 6- 



Note that in our notation R is always positive. For converging solutions, this implies 
that R is negative, growing towards zero. In our notation, therefore, converging solutions 
have negative (. This choice avoids the need to keep a ± sign in the expression of R for 
imploding and exploding solutions. 

The self similar solution exhibits time dependent characteristic pressure and density, 
which we denote as P(t) and N(t). For consistency with the BM notation we define the self 
similar variables as 

1 



7 2 M) = f 2 (t)g((), 
p(r,t) = P(t)/(C), 
n(r,t) = N(t)h(C)/9 1/2 (C) 



(5) 
(6) 
(7) 



where g, f and h are functions of the self similar variable ( and respectively describe the 
spatial profiles of the Lorentz factor, pressure and density. 



In terms of the self similar variables we have 



dt dT 



and 



24 + (r 2 -i/2)^ 



d_ = _ r2 ]_d_ 

dr Rd(' 



The energy conservation equation now becomes: 



(8) 



(9) 



= 2 (-2m + a - k)gf - (1 + 2(m + 1)0^/)' + /' 



(10) 



As in the work of BM, we are motivated by this equation to define 



X = l + 2(m + l)C, 



;n) 



to obtain 



= 2(-2m + a - k)gf - 2(m + 1) 



x-j-igf) - ~j-f 

dx dx . 



(12) 
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Here we have used the notation ^ = —m/2, ^ = — m — k, and ^ = — m/2 — k. The 
equation above shows that m and k must be constant to allow for a self similar flow. The 
relation between the derivatives of N, P, and T follows from the boundary conditions of 
strong relativistic shocks. Since we set x — 1 & t the shock, those boundary conditions also 
imply g(l) = /(l) = h(l) = 1 (see BM). 

Repeating the same procedure for the momentum equation will again result in equation 
(12). This is because the energy and momentum equations are identical to lowest order in 
1/r and our expansion discards any higher order terms. For this reason BM kept higher 
order terms. Instead, we find it simpler to use the difference equation between the energy 
and momentum conservation equations: 

+ ^ + ^"(4-^-2^ = 0, (13) 

Keeping only lowest order terms is again sufficient here. With the self similar variables 
equation 13 reads: 



= ( m + 1) 



- X g 2 — + 4(— g - —f) 
dx dx dx , 



+ (-m- k + 2a)g 2 f. (14) 



Similarly, for the particle conservation equation we have: 

^ = (-kg 2 - mg 2 + ag 2 - 2{m + lA/( X9 2 ~ 2g)/(m + 1). (15) 
dx dx 

We now solve equations (12), (14), and (15) as a set of linear equations to get df /dx, 
dg/dx, and dh/dx and obtain: 

1 dlogg (7m + 3k — 2a) — (m + a)xg 



gxdlogx (m+l)(g 2 x 2 -8gx + 4:) 



(16) 



1 dlogf _ 4(2m - q + fc)-(m + fc - 2a)gx ^ 



QXdlogx (m+l)(g 2 x 2 -Sgx + 4:) 

1 dlogh 2(9m + 5k - 4a) - 2(5m + 4k - 3a)gx + (m + k - a)g 2 x 2 
gxdlogx (m+ 1)(2 -gx)(g 2 X 2 ~ 8^x + 4) 

For the case a = 2 these reduce to the equations given in BM. 
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3. First type solutions 

The parameter m in the equations above must be found in an independent way. As is 
well known in the non-relativistic case, there are two types of self-similar solution. In the 
first type, m is found by global conservation laws. For arbitrary geometry, a, in the ultra 
relativistic regime this requirement reads 

E ~ 7 2 i?- fe i? 1+a ~ R -m-k+l+a _ ^ 

Therefore first type solutions have 

m — l + a — k. (20) 

Substituting this into the self similar equations and using the boundary conditions 
g(l) = f(l) = h(l) = 1 we obtain the very simple solution: 

9 = X' 1 (21) 

4fc-7-5a 

/ = x 3(2+ - fc) (22) 

and 

2fc-3-2c* 

h = x 2+a - k ■ (23) 

This is a generalized form of the BM solution, and reduces to the BM solution in spherical 
geometry where a — 2. 

First type solutions can be valid only if they contain a finite amount of energy. The 
energy in the solution is proportional to the integral 



/ 



4fc-7-5a 

fgdx~X 3(2+a - k) - (24) 



This is finite if 1 < x and < or if < x < 1 and fg^y > 0. 

1 < x if R(jn + 1) = R(2 + a — k) is positive. So the combination of the two possibilities 
above requires (4k — 7 — 5a)R < 0. Therefore, explosions diverging to infinite distances (R 
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is positive) can be of first type if k < (7 + 5a)/4. Converging solutions (where R is negative 
approaching zero) are of first type only for k > (7 + 5a) /A. 

For x < 1 these solutions are hollow and end at x = 0. The Lorentz factor diverges 
at x — 0) an d it takes infinite time for a fluid element to arrive from the shock (x = 1) to 
X = 0. Hollow solutions exist for diverging shocks in spherical symmetry with 4 < k < 17/4 
and for converging shocks in planar symmetry with 2 > k > 7/4. 

4. Second Type Solutions 

Second type solutions do not obey global conservation law. The true problem, therefore, 
cannot be completely described by a second type self similar solution. Those describe only 
part of the flow, in some region of interest while other regions deviate from the solutions. 
In order not to influence the self similar part, the non self similar parts must be separated 
from it by a sonic point, where the equations are singular. This requierment replaces the 
energy conservation as means of deducing the scaling of Lorentz factor with radius, i.e. 
finding m (see (5) for discussion of the non relativistic case and (7) for the relativistic case). 

To find m in second type solutions, we notice that the denominators of equations 16, 
17, and, 18 is independent of a, therefore the sonic line has the same value when expressed 
in terms of gx- The dependence of m on k can be easily found, using the same method as 
in Best and Sari. Looking for roots of the denominator of our equations we find the sonic 
line to be 

g X = 4 - 2V3. (25) 
While the fluid equations have other singular points 2 , physical considerations (see ref. (7)) 

2 The other singular points are gx = 2, which represents a fluid element maintaining a 
fixed x, and gx = 4 + 2^/3, which corresponds to the negative characteristic having a fixed 
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show that only the one given by (25) is the one we are looking for. This point separates 
the characteristic heading to the positive direction from the shock, which also goes in the 
positive direction per our definitions. Substituting this in the numerator of, for example, 
the equation of g and demanding that it vanish, we obtain: 

m = -2a(5 - 3"\/3) + (3 - 2y/l)k. (26) 

This reproduces the results of Best and Sari for a = 2. The hydrodynamic profiles can be 
found by substituting this value of m in the equations for g, f and h. The equations become 
very simple for this value of m since by definition the denominator and the numerator have 
a common factor. 



g = K 



1 — a 



-g X -2(2 + VS) 



(3-2v / 3)(fc-3a)/(a-l) 



-10a + 6aV3 + 3k - 2k^/3 + 1' 
where K is an arbitrary constant. The boundary conditions g(l) = 1 implies 

(3-2a/3) (fe-3a)/(a-l) 



- 1) + 4a( v / 3 - 1) - 2ky/3 + 4 + 2^3 



(a - 1) + 4a( v / 3 - 1) - 2k^3 + 4 + 2\/3 



(27) 



(28) 



For the pressure, we get 

U-2y^){k-M) 

log / = K + log 

a — 1 



- 1) - Aa + AaV3 - 2kV3 + A + 2V3] (29) 



or, with the boundary conditions /(l) = 1, 



./ 



^X(" - 1) - 4a + AaV3 - 2k^/3 + 4 + 2^3 
a - 1 - 4a + 4a v / 3 - 2^^ + 4 + 2^ 



(4-2 v / 3)(fe-3a)/(a-l) 



(30) 



For the density we get: 



log h = K- 



X- All those values are independent of the geometry, i.e. independent of the parameter a. 
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(2VS - 3) (2k - 1 + aV3 - 3a) (k - 3a) , r, - 
(-l-2a^ + W5-^ + a(l-a) '° S I' 1 " ° )9 * + * " ^ + ^ " ' ~ ^ 



-1 - 2a v / 3 + kV3 - V3 + a 
With the boundary conditions h — 1 at x — 1 we have: 



2a — k — ay3 , „ s , 

log( 9X -2). (31) 



(1 - + 4a - + - 4 - 2^3 

(1 - a) + 4a - 4a v / 3 + 2jb-s/3 - 4 - 2^ 



(2 v / 3-3)(2fc-l+QV / 3-3a)(fc-3c t ) 

( a -l-2 a y3 +fe V3-V3)( Q -l) -2 Q+fc + a ^ 

[2 — <7%] 0-1-20^+*:^-^ 

(32) 



For cylindrical geometry, a = 1, the above form is invalid and the solution (with the 
boundary condition of(?=latx = l) reads 

l0M = WTW> (33) 

log/= 3 ^|^(9X-l) (34) 

The condition we used to demand a smooth transition through the sonic point, which 
led to equation 26, is necessary but not sufficient. We still need to verify that the solution 
passes through the sonic point in the relevant range of the independent variable x- We 
therefore substitute gx = 4 — 2v^3 in the expression for g and find the value Xsonic- This 
is displayed in the figure below for the three possible values of a. We demand that Xsonic 
is within the range of values that x takes, i.e., we require Xsonic > 1 for diverging solution 
with m > — 1 or converging solutions with m < —1, and we require Xsonic < 1 fo r diverging 
solutions with m > — 1 or converging solutions with m < — 1. The results regarding the 
validity range of type-II solutions are given in table I and shown schematically in figure 1. 

A necessary condition for a second type solution to be correct, is that its energy 
decrease with time, that is, m > 1 + a — k for diverging solutions and the opposite for 
converging. We demonstrate this now for diverging solutions in spherical geometry: 

m u = -2a(5 - 3^) + (3 - 2V3)k >l + a-k = m I (35) 
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symmetry 


diverging 
type-I 


solutions (it — > oo) 
Type-II 


converging s 
Type-I 


olutions (R — > 0) 
Type-II 


planar 


k < 7/4 


fc > 2 


fc > 7/4 


fc < 1 


cylindrical 


k < 3 


fc > 3 


fc > 3 


fc < 3 


spherical 


k < 17/4 


A;>5 - a/3/4 ^4.13 


k > 17/4 


k < 4 



Table 1: The validity range of second type solutions. The choice of m according to equation 
26 guarantees that if a solution passes through a sonic point, the solution is smooth there. 
The conditions listed in the table verify that the solution does in fact pass though the sonic 
line within the physically relevant range of x- 

or 

k> (l + lla-6av / 3)/(4-2v / 3) (36) 

This gives k > 1 + a/3/4 for planar symmetry, k > 3 for cylindrical symmetry and 
k > 5 — a/3/4 for the spherical case. This is completely consistent with the above 
requirements. Indeed the requirement of decreasing energy follows from the condition that 
the solution passes through the sonic point: if the flow behind the sonic point cannot 
influence the flow ahead of the sonic point, energy could not flow from behind the sonic. 

5. The Type-I to Type-II Transition 

Note, that when mi = ran the two solutions listed above are the same, despite 
their seemingly different expressions. Clearly this must be the case since for a given m 
the solution to the self similar equations with the shock boundary conditions is unique. 
Nevertheless, it is instructive to verify this from the two explicit expressions given in the 
previous sections (e.g. equation 21 and 28). Let us consider the spherical case where 
mi = mn at k = 5 — a/3/4, and start by examining the second type solution. For 




Fig. 1. — Schematic description of the possible solutions for ID strong relativistic explosions. The 
three panels are for planar (top), cylindrical (middle) and spherical geometries (bottom). The x 
axis is the density index k (p oc r~ k ). Grey stripes indicate the region where first type solutions exist 
(k < 4 in the spherical case), and black stripes show where a second type solution exists. Hollow 
first type solutions, in which all the shocked fluid is confined between the shock and some internal 
surface are indicated by a hollow grey stripe. The y-axis of grey and black stripes is arbitrary. The 
solid line is the position of the sonic point, Xsonic- 



-14- 



k m 5 — a/3/4, the numerator of equation (28) is close to zero. This means than gx has 
to be close to unity to make the numerator close to zero as well. This agrees with the 
expression for g in the first type solution g = x 1 - 

Where is the sonic point in type-II solutions for k close to k = 5 — a/3/4? Expanding 
equation (28) we obtain 

' 1 (37) 



_t-(5- v^/4)_ 

So that (? SO mc diverges when k is close to fc = 5 — a/3/4. Since at the sonic point 
dsonicXsonic = 4 — 2 a/3 is finite, Xsonic approaches zero. This is very similar to the behavior of 
Type-I solutions as those end at x — with divergence of g. Nevertheless, Type-II solutions 
continue past x — into negative x an d they are not hollow. 

An unfortunate property of the solutions around k — 4 (in spherical geometry) is 
that the profile approaches non relativistic temperatures very quickly as a function of the 
distance from the shock. We define non relativistic temperatures to occur where p/n < 1. 
We have 

| ~ T^- ~ r(l + 2(4 - k)C)$=V ~ 7exp (-8C/3) (38) 



So the temperature becomes non relativistic at 

Gvi^lnr. (39) 

o 

This is just a few times R/T 2 behind the shock, even for quite large T. For k ~ 4 the 
solution extends many times -R/7 2 behind the shock. We therefore expect that for moderate 
values of T, the ultra relativistic self similar solution gives an approximate rather than an 
accurate description of the flow. 

Clearly it is of importance to understand which of the two solutions apply. In principle, 
one can attempt this with numerical simulations. However, the problem of non relativistic 
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temperatures noted above makes numerical investigation difficult as it requires extremely 
large Lorentz factors which are hard to achieve. 

6. Discussion 

We have explored the possible ultra-relativistic self similar solutions in planar, 
cylindrical, and spherical geometries, containing diverging or converging shocks, and 
allowing for a general powerlaw density profile p oc r~ k . Relativistic implosions in all three 
geometries and with a general powerlaw density profile are treated here for the first time. 
While this paper was being prepared, an independent paper appeared treating imploding 
relativistic shock waves for spherical geometry with constant density (10). 

We show that as a function of the density powerlaw index k, there is a sequence of self 
similar solutions. In the explosive scenarios, the sequence is from first type solutions to 
hollow first type solutions and then second type solutions as k increases. The sequence is 
reversed for imploding solutions. Those are of second type for low k and first type for large 
k. In both cases, explosions and implosions, the solution is of second type if the density 
ahead of the shock decreases sufficiently fast and of first type if it increases sufficiently 
fast. The physical interpretation of this is clear. A shock wave that propagates into a 
density profile that decreases sufficiently fast will accelerates and "run away" from the 
fluid behind it. The shock becomes "causally disconnected" from the downstream fluid, 
i.e. a second type solution. On the other hand, if the shock wave does not accelerate 
enough, sound waves from behind it will be able to catch up and the whole flow is causally 
connected, i.e. a first type solution. In some of the geometries (converging planar case and 
diverging spherical case) first type solutions become hollow before they turn into second 
type solutions. This sequence of self similar solutions and its physical origin apply to the 
non relativistic case as well. 
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Our analysis raises some interesting riddles. In some cases (planar diverging shocks or 
spherical converging shocks) this sequence has a "gap" between the first and second type 
solutions. In other cases (planar converging shocks and spherical diverging shocks) there 
seems to be an overlap in the sequence between first and second type solutions: current 
considerations allow for the existence of both self similar solutions for the same value of k. 

Recently, Gruzinov (11) suggested a solution in the "gap" region in the non relativistic 
case which he called a "third type" self similar solution. In his third type solution, the 
infinite mass located at the origin acts like a piston moving at constant velocity. Preliminary 
investigation seems to support a relativistic analog of Gruzinov's analysis. In the relativistic 
case, however, the solutions ara more interesting since the piston with large inertia does 
not move at a constant speed but accelerates. Further research is needed to establish this 
connection. 

This research was partially supported by a NASA ATP grant. RS is a Packard Fellow 
and an Alfred P. Sloan Research Fellow. We thank Margaret Pan and Chris Matzner for 
helpful discussions. 
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Fig. 2. — First and second type solutions for k = 5— -y/3/4+0.01 plotted as a function ofC- The top 
plot shows the Lorentz factor squared profile g and the bottom plot the pressure profile f. The two 
solutions are very similar close to the shock (Q = 0). The Type-I solution is hollow, so the pressure 
drops to zero and 7 diverges at some finite distance from the shock Q = — l/(2m + 2) = 2 + 
The second type solution continues smoothly but does show a relatively sharp drop around the same 
distance. 
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; 

Fig. 3. — Comparison of numerical integration of the time dependent hydrodynamical equations 
and the prediction of the self similar solution. Plotted, top to bottom, are the Lorentz factor, the 
pressure and density relative to that just behind the shock, for an imploding spherical shockwave 
into matter of constant density: a = 2, k = 0. 
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